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Abstract 

As first discovered by Choptuik, the black hole threshold in the space of initial data for gen- 
eral relativity shows both surprising structure and surprising simplicity. Universality, power- 
law scaling of the black hole mass, and scale echoing have given rise to the term "critical 
phenomena". They are explained by the existence of exact solutions which are attractors 
within the black hole threshold, that is, attractors of codimension one in phase space, and 
which are typically self-similar. Critical phenomena give a natural route from smooth initial 
data to arbitrarily large curvatures visible from infinity, and are therefore likely to be relevant 
for cosmic censorship, quantum gravity, astrophysics, and our general understanding of the 
dynamics of general relativity. 
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1 Introduction 

1.1 Overview of the subject 

Take generic initial data in general relativity, adjust any one parameter p of the initial data to 
the threshold of black hole formation, and compare the resulting spacetimes as a function of p. In 
many situations, the following critical 'phenomena are then observed: 

• Near the threshold, black holes with arbitrarily small masses can be created, and the black 
hole mass scales as 

M«(p-p,p, (1) 
where p parameterises the initial data and black holes form for p > . 

• The critical exponent 7 is universal with respect to initial data, that is, independent of the 
particular 1-parameter family, although it depends on the type of collapsing matter. 

• In the region of large curvature before black hole formation, the spacetime approaches a 
self-similar solution which is also universal with respect to initial data, the critical solution. 
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Critical phenomena were discovered by Choptuik |46j in numerical simulations of a spherical 
scalar field. They have been found in numerous other numerical and analytical studies in spherical 
symmetry, and a few in axisymmetry; in particular critical phenomena have been seen in the 
collapse of axisymmetric gravitational waves in vacuum [J. Scaling laws similar to the one for 
black hole mass have been discovered for black hole charge and conjectured for black hole angular 
momentum. 

It is still unclear how universal critical phenomena in collapse are with respect to matter types 
and beyond spherical symmetry, in particular for vacuum collapse. Progress is likely to be made 
over the next few years with better numerical simulations. 

Critical phenomena can usefully be described in dynamical systems terms. A critical solution is 
then characterised as an attracting fixed point within a surface that divides two basins of attraction, 
a critical surface in phase space. Such a fixed point can be either a stationary spacetime, or one 
that is scale-invariant and self-similar. The latter is relevant for the ("type 11") critical phenomena 
sketched above. We shall also see how the dynamical systems approach establishes a connection 
with critical phenomena in statistical mechanics. 

Therefore, we could define the field of critical phenomena in gravitational collapse as the study 
of the boundaries among the basins of attraction of different end states of self-gravitating systems, 
such as black hole formation or dispersion. In our view, the main physical motivation for this 
study is that those critical solutions which are self-similar provide a way of achieving arbitrarily 
large spacetime curvature outside a black hole, and in the limit a naked singularity, by fine-tuning 
generic initial data for generic matter to the black hole threshold. Those solutions are therefore 
likely to be important for quantum gravity and cosmic censorship. 

1.2 Plan of this review 

Faced with more material to review, we have attempted to make this update shorter and more 
systematic than the original 1999 version. We begin with the abstract theory in Sec. [21 Since 
1999, progress on the theory side has mainly been made on the global spacetime structure of the 
critical solution and cosmic censorship in the spherical scalar field model. We have included this 
new material in an enlarged Sec. [3] on the spherical scalar field, although we hope it will turn out 
to be sufficiently generic to merit inclusion in Sec. [21 Nonspherical perturbations of the spherical 
scalar field are also discussed in Sec. [31 

In Sec. [H we review the rich phenomenology that has been found in many other systems re- 
stricted to spherical symmetry. Numerical work in spherical symmetry has proliferated since 1999, 
but we have tried to keep this section as short as possible. There has been less progress in going 
beyond spherical symmetry than we anticipated in 1999, even though we continue to believe that 
important results await there. What is known today is summarised in Sec. [5l 

The reader unfamiliar with the topic is advised to begin with either Sees. 12.1112.31 which give 
the key theory of universality, self-similarity and scaling, or Secs. l3.lll3.2i which describe the classic 
example, the massless scalar field. 

This review is limited to numerical and theoretical work on phenomena at the threshold of black 
hole formation in 3-|-l-dimensional general relativity. We report only briefly on work in higher 
and lower spacetime dimensions and non-gravity systems that may be relevant as toy models for 
general relativity. We exclude other work on self-similarity in general relativity and work on critical 
phenomena in other areas of physics. 

Other reviews on the subject are [l24], [13], [92], [96], [48], [49], [32], [98], [HE]- The 2002 
review [lOlj by Gundlach gives more detailed explanations on some of the basic aspects of the 
theory. 
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2 Theory 



In this Section we describe the basic theory underlying critical collapse of the type that forms 
arbitrarily small black holes (later called type II, see also Sec. l2.4p . We begin with the mathematical 
origin of its three main characteristics, which were already summarised in the introduction: 

• universality with respect to initial data; 

• scale-invariance of the critical solution; 

• black hole mass scaling. 

2.1 Universality 

Consider GR as an infinite-dimensional continuous dynamical system. Points in the phase space 
are initial data sets (3-metric, extrinsic curvature and suitable matter variables, which together 
obey the Einstein constraints). We evolve with the Einstein equations in a suitable gauge (see 
Sec. 12. Sp . Solution curves of the dynamical system are spacetimes obeying the Einstein- matter 
equations, sliced by specific Cauchy surfaces of constant time t. 

An isolated system in GR can end up in qualitatively different stable end states. Two possibil- 
ities are the formation of a single black hole in collapse, or complete dispersion of the mass-energy 
to infinity. For a massless scalar field in spherical symmetry, these are the only possible end states 
(see Sec. |3]). Any point in phase space can be classified as ending up in one or the other type of 
end state. The entire phase space therefore splits into two halves, separated by a "critical surface" . 

A phase space trajectory that starts on a critical surface by definition never leaves it. A critical 
surface is therefore a dynamical system in its own right, with one dimension fewer than the full 
system. If it has an attracting fixed point, such a point is called a critical point. It is an attractor 
of codimension one in the full system, and the critical surface is its attracting manifold. The fact 
that the critical solution is an attractor of codimension one is visible in its linear perturbations: it 
has an infinite number of decaying perturbation modes spanning the tangent plane to the critical 
surface, and a single growing mode not tangential to it. 

As illustrated in Figs. [T] and [2l any trajectory beginning near the critical surface, but not 
necessarily near the critical point, moves almost parallel to the critical surface towards the critical 
point. Near the critical point the evolution slows down, and eventually moves away from the 
critical point in the direction of the growing mode. This is the origin of universality. All details 
of the initial data have been forgotten, except for the distance from the black hole threshold. The 
closer the initial phase point is to the critical surface, the more the solution curve approaches the 
critical point, and the longer it will remain close to it. We should stress that this phase picture 
is extremely simplified. Some of the problems associated with this simplification are discussed in 
Sec. [231 

2.2 Self-similarity 

Fixed points of dynamical systems often have additional symmetries. In the case of type II critical 
phenomena, the critical point is a spacetime that is self-similar, or scale-invariant. These sym- 
metries can be discrete or continuous. The critical solution of a spherically symmetric perfect 
fluid (see Sec. 14. 2|) . has continuous self-similarity (CSS). A CSS spacetime is one that admits a 
homothetic vector field defined by [36] : 





(2) 



4 




Figure 1: The phase space picture for the black hole threshold in the presence of a critical point. 
Every point correspond to an initial data set, that is, a 3-metric, extrinsic curvature, and matter 
fields. (In type II critical collapse these are only up to scale). The arrow lines are solution curves, 
corresponding to spacetimes, but the critical solution, which is stationary (type I) or self-similar 
(type II) is represented by a point. The line without an arrow is not a time evolution, but a 1- 
parameter family of initial data that crosses the black hole threshold at p = p*. The 2-dimensional 
plane represents an (oo — l)-dimensional hypersurface, but the third dimension represents really 
only one dimension. 
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Figure 2: A diflFerent phase space picture, specifically for type II critical collapse, and two 2- 
dimensional projections of the same picture. In contrast with Fig. [U one dimension of the two 
representing the (infinitely many) decaying modes has been suppressed. The additional axis now 
represents a global scale which was suppressed in Fig. [U so that the scale-invariant critical solution 
CS is now represented as a straight line (in red). Several members of a family of initial conditions 
(in blue) are attracted by the critical solution and then depart from it towards black hole formation 
(A or B) or dispersion (D). Perfectly fine-tuned initial data asymptote to the critical solution with 
decreasing scale. Initial conditions starting closer to perfect fine tuning produce smaller black holes, 
such that the parameter along the line of black hole endstates is —IuMbh- Two 2-dimensional 
projections of the same picture are also given. The horizontal projection of this figure is the same 
as the vertical projection of Fig. [TJ 



6 



In coordinates x'^ = {t,x^) adapted to the symmetry, so that 



(3) 



the metric coefficients are of the form 

.gA<i^(T, x') = e^'^^gf,^{x'), (4) 

where the coordinate r is the negative logarithm of a spacetime scale, and the remaining three 
coordinates x^ can be thought of angles around the singular spacetime point t = oo (see Sec. 13. 3p . 

The critical solution of other systems, in particular the spherical scalar field (see Sec. [3]) and 
axisymmetric gravitational waves (see Sec. 15. 2p . show discrete self-similarity (DSS). The simplest 
way of defining DSS is in adapted coordinates, where 

9t.u{T, x') = e"2^g^^(r, x'), (5) 

such that g^i,j{T, x') is periodic in t with period A. More formally, DSS can be defined as a discrete 
conformal isometry [91j . 

Using the gauge freedom of general relativity, the lapse and shift in the ADM formalism can 
be chosen (non-uniquely) so that the coordinates become adapted coordinates if and when the 
solution becomes self-similar (see Sec. 12. 5|) . t is then both a time coordinate (in the usual sense 
that surfaces of constant time are Cauchy surfaces), and the logarithm of overall scale at constant 
x^ . The minus sign in ^ and hence ([U and (O, is a convention assuming that smaller scales are 
in the future. The time parameter used in Fig. [2] is of this type. 



2.3 Mass scaling 

Let Z stand for a set of scale-invariant variables of the problem, such as cjf^i, and suitably rescaled 
matter variables. If the dynamics is scale-invariant (this is the case exactly for example for the 
scalar field, and approximately for other systems, see Sec. 12. 6|) . then Z{x) is an element of the 
phase space factored by overall scale, and Z{x, r) a solution. Note that Z{x) is an initial data set 
for GR only up to scale. The overall scale is supplied by r. 

For simplicity, assume that the critical solution is CSS. It can then be written as Z{x,t) — 
Zf,{x). Its linear perturbations can depend on t only exponentially. To linear order, the solution 
near the critical point must be of the form 



Z(a;,T) ~Z*(a;)+^Q(p)e^-"Z,(a;). (6) 



The perturbation amplitudes Ci depend on the initial data, and hence on p. As is a critical 
solution, by definition there is exactly one A; with positive real part (in fact it is purely real), say 
Aq. As r — > oo, all other perturbations vanish. In the following we consider this limit, and retain 
only the one growing perturbation. 

From our phase space picture, the evolution ends at the critical solution for p ~ p,,, so we must 
have Co(p*) — 0. Linearising in p around p*, we obtain 

lim Z(x,T)~Z,(x) + ^(p-p,)e^«^Zo(x). (7) 

T^OO dp 

For p 7^ p*, but close to it, the solution has the approximate form ([7]) over a range of r. Now 
we extract Cauchy data at one particular p-dependent value of t within that range, namely 
defined by 

^5^(p-p.)e-^-^., (8) 
dp 
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where e is some constant ^ 1 such that at this r the hnear approximation is stiU vahd. At 
sufficiently large r, the linear perturbation has grown so much that the linear approximation 
breaks down, and for Co > a black hole forms while for Cq < the solution disperses. The 
crucial point is that we need not follow this evolution in detail, nor does the precise value of e 
matter. It is sufficient to note that the Cauchy data at t = are 

Z{x,T^) ~ Z^{x) + eZo{x). (9) 

Due to the funnelling effect of the critical solution, the data at is always the same, except for an 
overall scale, which is given by e~'^' . For example, the physical spacetime metric, with dimension 
(length)^ is given by g^^, = e~'^'^g^^, and similar scalings hold for the matter variables according 
to their dimension. In particular, as e""^* is the only scale in the initial data ([5]), the mass of the 
final black hole must be proportional to that scale. Therefore 

M cx e"^* cx (10) 

and, comparing with Eq. we have found the critical exponent 7 = I/Aq. 

When the critical solution is DSS, a periodic or fine structure of small amplitude is superimposed 
on this basic power law [94l 1122] : 

In i\/ = 7 ln(p - p*) + c + / (7 ln(p - p,) + c) , (11) 

where f{z) has period A and is universal, and only c depends on the initial data. As the critical 
solution is periodic in r with period A the number N of scaling "echos" is approximated by 

TV ~ A"^7ln|p-p*| +const. (12) 

Note that this holds for both supercritical and subcritical solutions. 

2.4 Type I 

In type I critical phenomena, the same phase space picture as in Sec. 12.11 applies, but the critical 
solution is now stationary or time-periodic instead of self-similar or scale-periodic. It also has 
a finite mass and can be thought of as a metastable star. (Type I and II were so named after 
first and second order phase transitions in statistical mechanics, in which the order parameter is 
discontinuous and continuous, respectively.) Universality in this context implies that the black 
hole mass near the threshold is independent of the initial data, namely a certain fraction of the 
mass of the stationary critical solution. The dimensionful quantity that scales is not the black 
hole mass, but the lifetime tp of the intermediate state where the solution is approximated by the 
critical solution. This is clearly 

tp = — In |p — p*| + const. (13) 
Ao 

Type I critical phenomena occur when a mass scale in the field equations becomes dynamically 
relevant. (This scale does not necessarily set the mass of the critical solution absolutely: there 
could be a family of critical solutions selected by the initial conditions.) Conversely, as the type II 
power law is scale-invariant, type II phenomena occur in situations where either the field equations 
do not contain a scale, or this scale is dynamically irrelevant. Many systems, such as the massive 
scalar field, show both type I and type II critical phenomena, in different regions of the space of 
initial data [33] . 
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2.5 Coordinate choices for the dynamical systems picture 

The time evolution of Cauchy data in GR can only be considered as a dynamical system if the 
ADM evolution equations are complemented by a prescription for the lapse and shift. To realise the 
phase space picture of Sec. 12.11 the critical solution must be a fixed point or limit cycle. We have 
seen how coordinates adapted to the self-similarity can be constructed, but is there a prescription 
of the lapse and shift for arbitrary initial data, such that, given initial data for the critical solution, 
the resulting time evolution actively drives the metric to a form ^ that explicitly displays the 
self-similarity? 

Garfinkle and Gundlach [82] have suggested several combinations of lapse and shift conditions 
that leave CSS spacetimes invariant and turn the Choptuik DSS spacetime into a limit cycle. (See 
[871 178j for partial successes.) Among these, the combination of maximal slicing with minimal 
strain shift has been suggested in a different context but for related reasons [176) . Maximal slicing 
requires the initial data slice to be maximal (/^a" — 0), but other prescriptions, such as freezing 
the trace of K together with minimal distortion, allow for an arbitrary initial slice with arbitrary 
spatial coordinates. 

All these coordinate conditions are elliptic equations that require boundary conditions, and 
will turn CSS spacetimes into fixed points (or DSS into limit cycles) only given correct boundary 
conditions. Roughly speaking, these boundary conditions require a guess of how far the slice is 
from the accumulation point t ~ t^, and answers to this problem only exist in spherical symmetry. 
Appropriate boundary conditions are also needed if the dynamical system is extended to include 
the lapse and shift as evolved variables, turning the elliptic equations for the lapse and shift into 
hyperbolic or parabolic equations. 

Turning a CSS or stationary spacetime into a fixed point of the dynamical system also requires 
an appropriate choice of the phase space variables Z{x'^). To capture CSS (or DSS) solutions, one 
needs scale-invariant variables. Essentially, these can be constructed by dimensional analysis. The 
coordinates a;* and r are dimensionless, Ze~^ has dimension length, and has dimension P. The 
scaling for the ADM and any matter variables follows. 

Even with a prescription for the lapse and shift in place, a given spacetime does not correspond 
to a unique trajectory in phase space. Rather, for each initial slice through the same spacetime one 
obtains a different slicing of the entire spacetime. A possibility for avoiding this ambiguity would 
be to restrict the phase space further, for example by restricting possible data sets to maximal or 
constant extrinsic curvature slices. 

Another open problem is that in order to talk about attractors and repellcrs on the phase space 
we need to define a norm on a suitable function space which includes both asymptotically flat data 
and data for the exact critical solution. The norm itself must favour the central region and ignore 
what is further out and asymptotically flat if all black holes of the same mass are to be considered 
as the same endstate. 

2.6 Approximate self-similarity and universality classes 

The field equations for the massless scalar field coupled to the Einstein equations are scale-free. 
Realistic matter models introduce length scales, and the field equations then do not allow for exactly 
self-similar solutions. They may however admit solutions which are CSS or DSS asymptotically on 
small spacetime scales as the dimensionful parameters become irrelevant, including type II critical 
solutions [47l [33l [50] . This can be explored by a formal expansion in powers of the small parameter 
Le~'^ , where L is a parameter with dimensions length in the evolution equations. The zeroth order 
of the expansion is the self-similar critical solution of the system with L = 0. A similar ansatz 
can be made for the linear perturbations of the resulting background. The zeroth order of the 
background expansion determines A exactly and independently of L, and the zeroth order term 
of the linear perturbation expansion determines the critical exponent I/Aq exactly, so that there 
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is no need in practice to calculate any higher orders in L to make predictions for type II critical 
phenomena where they occur. (With L 7^ 0, the basin of attraction of the type II critical solution 
will depend on L, and type I critical phenomena may also occur, see Sec. 12.41 A priori, there could 
also be more than one type II critical solution for L = 0, although this has not been observed.) 

This procedure has been carried out for the Einstein- Yang-Mills system [93] and for massless 
scalar electrodynamics |102j . Both systems have a single length scale 1/e (in geometric units 
c = G = 1), where e is the gauge coupling constant. All values of e can be said to form one 
universality class of field equations [108] represented by e = 0. This notion of universality classes 
is fundamentally the same as in statistical mechanics. Other examples include modifications to 
the perfect fluid equation of state that do not affect the limit of high density [152] . A simple 
example is that any scalar field potential V{(t>) becomes dynamically irrelevant compared to the 
kinetic energy |V(/)p in a self-similar solution [47], so that all scalar fields with potentials are in 
the universality class of the free massless scalar field. Surprisingly, even two different models like 
the SU{2) Yang-Mills and SU{2) Skyrme models in spherical symmetry are members of the same 
universality class |20). 

If there are several scales Lq, Li, L2 etc. present in the problem, a possible approach is to set 
the arbitrary scale in (|28p equal to one of them, say Lq, and define the dimensionless constants 
li = Li/Lo from the others. The scope of the universality classes depends on where the U appear in 
the field equations. If a particular Li appears in the field equations only in positive integer powers, 
the corresponding li appears only multiplied by e~'^, and will be irrelevant in the scaling limit. All 
values of this U therefore belong to the same universality class. From the example above, adding a 
quartic self- interaction X(j)* to the massive scalar field gives rise to the dimensionless number X/m? 
but its value is an irrelevant (in the language of renormalisation group theory) parameter. 

Contrary to the statement in [102 , we conjecture that massive scalar electrodynamics, for any 
values of e and m, is in the universality class of the massless uncharged scalar field in a region of 
phase space where type II critical phenomena occur. Examples of dimensionless parameters which 
do change the universality class are the k of the perfect fluid, the k of the 2-dimensional sigma 
model or, probably, a conformal coupling of the scalar field [45j (the numerical evidence is weak 
but a dependence should be expected). 

2.7 The analogy with critical phase transitions 

Some basic aspects of critical phenomena in gravitational collapse, such as fine-tuning, universality, 
scale-invariant physics, and critical exponents for dimensionful quantities, can also be identified in 
critical phase transitions in statistical mechanics (see [193] for an introductory textbook). 

From an abstract point of view, the objective of statistical mechanics is to derive relations 
between macroscopic observables A of the system and macroscopic external forces / acting on it, 
by considering ensembles of microscopic states of the system. The expectation values (A) can be 
generated as partial derivatives of the partition function 



Here the fi are parameters of the Hamiltonian such as the strength of intermolecular forces, and / 
are macroscopic quantities which are being controlled, such as the temperature or magnetic field. 

Phase transitions in thermodynamics are thresholds in the space of external forces / at which 
the macroscopic observables A, or one of their derivatives, change discontinuously. We consider 
two examples: the liquid-gas transition in a fluid, and the ferromagnetic phase transition. 

The liquid-gas phase transition in a fluid occurs at the boiling curve p — Pb{T). In crossing 
this curve, the fluid density changes discontinuously. However, with increasing temperature, the 





(14) 



mic restates 



10 



difference between the liquid and gas density on the boihng curve decreases, and at the the critical 
point (p, =P6(T*),T*) it vanishes as a non-integer power: 

Pliquid - Pgas ~ - T)^. (15) 

At the critical point an otherwise clear fluid becomes opaque, due to density fluctuations appearing 
on all scales up to scales much larger than the underlying atomic scale, and including the wavelength 
of light. This indicates that the fluid near its critical point is approximately scale-invariant (for 
some range of scales between the size of molecules and the size of the container). 

In a ferromagnetic material at high temperatures, the magnetisation m of the material (align- 
ment of atomic spins) is determined by the external magnetic field B. At low temperatures, the 
material shows a spontaneous magnetisation even at zero external field. In the absence of an ex- 
ternal field this breaks rotational symmetry: the system makes a random choice of direction. With 
increasing temperature, the spontaneous magnetisation m decreases and vanishes at the Curie 
temperature as 

|m| ~ (T, - T)T. (16) 

Again, the correlation length, or length scale of a typical fluctuation, diverges at the critical point, 
indicating scale-invariant physics. 

Quantities such as |m| or puquid — Pgas are called order parameters. In statistical mechanics, 
one distinguishes between first-order phase transitions, where the order parameter changes discon- 
tinuously, and second-order, or critical, ones, where it goes to zero continuoiisly. One should think 
of a critical phase transition as the critical point where a line of first-order phase transitions ends 
as the order parameter vanishes. This is already clear in the fluid example. In the ferromagnet 
example, at first one seems to have only the one parameter T to adjust. But in the presence of 
a very weak external field, the spontaneous magnetisation aligns itself with the external field B, 
while its strength is to leading order independent of B. The function m(B, T) therefore changes 
discontinuously at B = 0. The line B = for T < is therefore a line of first order phase transi- 
tions between directions (if we consider one spatial dimension only, between m up and m down). 
This line ends at the critical point (B = 0, T = T,) where the order parameter |m| vanishes. The 
critical value B = of B is determined by symmetry; by contrast depends on microscopic 
properties of the material. 

We have already stated that a critical phase transition involves scale-invariant physics. In 
particular, the atomic scale, and any dimcnsionful parameters associated with that scale, must 
become irrelevant at the critical point. This is taken as the starting point for obtaining properties 
of the system at the critical point. 

One first defines a semi-group acting on micro-states: the renormalisation group. Its action 
is to group together a small number of adjacent particles as a single particle of a fictitious new 
system by using some averaging procedure. This can also be done in a more abstract way in Fourier 
space. One then defines a dual action of the renormalisation group on the space of Hamiltonians 
by demanding that the partition function is invariant under the renormalisation group action: 

E .E (17) 

microstatcs microstatcs' 

The renormalised Hamiltonian is in general more complicated than the original one, but it can be 
approximated by the same Hamiltonian with new values of the parameters ^, and external forces /. 
(At this stage it is common to drop the distinction between /i and /, as the new /i' and /' depend 
on both ix and /.) Fixed points of the renormalisation group correspond to Hamiltonians with the 
parameters at their critical values. The critical values of many of these parameters will be or 
oo, meaning that the dimensionful parameters they were originally associated with are irrelevant. 
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Because a fixed point of tlie renormalisation group can not have a preferred length scale, the only 
parameters that can have nontrivial values are dimensionless. 

The behaviour of thermodynamical quantities at the critical point is in general not trivial to 
calculate. But the action of the renormalisation group on length scales is given by its definition. 
The blowup of the correlation length ^ at the critical point is therefore the easiest critical exponent 
to calculate. The same is true for the black hole mass, which is just a length scale. We can 
immediately reinterpret the mathematics of Section 12.31 as a calculation of the critical exponent 
for ^, by substituting the correlation length ^ for the black hole mass M , — T for p — p*, and 
taking into account that the r-evolution in critical collapse is towards smaller scales, while the 
renormalisation group flow goes towards larger scales: f therefore diverges at the critical point, 
while M vanishes. 

In type II critical phenomena in gravitational collapse, we should think of the black hole mass 
as being controlled by the functions P and Q on phase space defined by Eq. (|26p. Clearly, P is 
the equivalent of the reduced temperature T — T^. Gundlach [lOOj has suggested that the angular 
momentum of the initial data can play the role of B, and the final black hole angular momentum 
the role of m. Like the magnetic field, angular momentum is a vector, with a critical value that 
must be zero because all other values break rotational symmetry. 

3 The scalar field 

Critical phenomena in gravitational collapse were first discovered by Choptuik [45t l46l |47] in the 
model of a spherically symmetric, massless scalar field minimally coupled to general relativity. 
The scalar field matter is both simple, and acts as a toy model in spherical symmetry for the effects 
of gravitational radiation. Given that it is still the best-studied model in spherical symmetry, we 
review it here as a case study. For other numerical work on this model, see |1Q5I [77l I107[ [74l 
I171J . Important analytical studies of gravitational collapse in this model have been carried out by 
Christodoulou [55l [Ml [SH [Ml IMl [60l E] . 

We first review the field equations and Choptuik's observations at the black hole threshold, 
mainly as a concrete example for the general ideas discussed above. We then summarise more 
recent work on the global structure of Choptuik's critical solution, which throws an interesting light 
on cosmic censorship. In particular, the exact critical solution contains a curvature singularity that 
is locally and globally naked, and any critical solution obtained in the limit of perfect fine-tuning of 
asymptotically flat initial data is at least locally naked. By perturbing around spherical symmetry, 
the stability of the Choptuik solution in the full phase space can be investigated, and the scaling 
of black hole angular momentum can be predicted. By embedding the real scalar field in scalar 
electrodynamics and perturbing around the Choptuik solution, the scaling of black hole charge can 
be predicted. 

3.1 Field equations in spherical symmetry 

The Einstein equations are 

Gat - 8^ (^V„0Vb0 - ig„bV,0V=0^ , (18) 

and the matter equation is 

V„V> = 0. (19) 

Note that the matter equation of motion is contained within the contracted Bianchi identities. 
Choptuik chose Schwarzschild-like coordinates 

ds^ = -a^ (r, t) df + (r, t) dr^ + dVl? , (20) 
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where dfi'^ = d9^ + sin^ 9 dip^ is the metric on the unit 2-sphere. This choice of coordinates is 
defined by the radius r giving the surface area of 2-spheres as 47rr^, and by t being orthogonal 
to r (polar-radial coordinates). One more condition is required to fix the coordinate completely. 
Choptuik chose a = 1 at r = 0, so that t is the proper time of the central observer. 
In the auxiliary variables 

* = n--0.t, (21) 

a ' 

the wave equation becomes a first-order system, 

^. = (-n) , 



r 



(22) 



In spherical symmetry there are four algebraically independent components of the Einstein equa- 
tions. Of these, one is a linear combination of derivatives of the other and can be disregarded. The 
other three contain only first derivatives of the metric, namely a_t, a^r and a^r, and are 



a^r a — 1 
a r a r — 1 



a a r 



= 27rr(n2 + $2)^ (23) 

= 0, (24) 
= 47rr$n. (25) 



Because of spherical symmetry, the only dynamics is in the scalar field equations ()22p . The metric 
can be found by integrating the ODEs (|23p and ([24|) for a and a at any fixed t, given (j) and 11. 
Eq. (P5|) can be ignored in this "fully constrained" evolution scheme. 



3.2 The black hole threshold 

The free data for the system are the two functions n(r, 0) and $(r, 0). Choptuik investigated 
several one-parameter families of such data by evolving the data for many different values of the 
parameter. Simple examples of such families are n(r, 0) = and a Gaussian for $(r, 0), with the 
parameter p taken to be either the amplitude of the Gaussian, with the width and centre fixed, 
or the width, with position and amplitude fixed, or the position, with width and amplitude fixed. 
For sufficiently small amplitude (or the peak sufficiently wide), the scalar field will disperse, and 
for sufficiently large amplitude it will form a black hole. 

Generic one-parameter families behave in this way, but this is difficult to prove in generality. 
Christodoulou showed for the spherically symmetric scalar field system that data sufficiently weak 
in a well-defined way evolve to a Minkowski-like spacetime [56j , and that a class of sufficiently 
strong data forms a black hole ^ . 

Choptuik found that in all 1-parameter families of initial data he investigated he could make 
arbitrarily small black holes by fine-tuning the parameter p close to the black hole threshold. An 
important fact is that there is nothing visibly special to the black hole threshold. One cannot tell 
that one given data set will form a black hole and another one infinitesimally close will not, short 
of evolving both for a sufficiently long time. 

As p —^ along the family, the spacetime varies on ever smaller scales. Choptuik developed 
numerical techniques that recursively refine the numerical grid in spacetime regions where details 
arise on scales too small to be resolved properly. In the end, he could determine up to a 
relative precision of 10~^^, and make black holes as small as 10~^ times the ADM mass of the 
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spacetime. The power-law scaling ([T0|) was obeyed from those smallest masses up to black hole 
masses of, for some families, 0.9 of the ADM mass, that is, over six orders of magnitude [37]. 
There were no families of initial data which did not show the universal critical solution and critical 
exponent. Choptuik therefore conjectured that 7 is the same for all one-parameter families of 
smooth, asymptotically flat initial data that depend smoothly on the parameter, and that the 
approximate scaling law holds ever better for arbitrarily small p — p^,. 

It is an empirical fact that typical one-parameter families cross the threshold only once, so 
that there is every indication that it is a smooth submanifold, as we assumed in the phase space 
picture. Taking into account the discussion of mass scaling above, we can formally write the black 
hole mass as a functional of the initial data z ~ {(t>{r, 0), n(r, 0)) exactly as 

M[z] = Q[z]H{P[z]){P[z]y, (26) 

where P and Q are smooth functions on phase space and H is the Heaviside function. (Q could be 
absorbed into P). 

In hindsight, polar-radial gauge is well-adapted to self-similarity. In this gauge, discrete self- 
similarity corresponds to 

Z{r, t)^Z (e"^r, e^H) (27) 

for any integer n, where Z stands for any one of the dimensionless quantities a, a or (and 
therefore also for rll or r$). With 

. = -^, . = -ln(-i^), t<t„ (28) 

discrete self-similarity is 

Z{x,T + I^)^ Z{x,t). (29) 

The dimensionful constants and L depend on the particular one-parameter family of solutions, 
but the dimensionless critical fields a*, a* and 0*, and in particular their dimensionless period 
A, are universal. Empirically, A ~ 3.44 for the scalar field in numerical time evolutions, and 
A = 3.445452402(3) from a numerical construction of the critical solution based on exact self- 
similarity and analyticity [149] . 

3.3 Global structure of the critical solution 

In adapted coordinates, the metric of the critical spacetime is of the form e^^"^ times a regular 
metric. From this general form alone, one can conclude that r = 00 is a curvature singularity, 
where Riemann and Ricci invariants blow up like e*"^ (unless the spacetime is flat), and which is 
at finite proper time from regular points in its past. The Weyl tensor with index position C^^cd 
is conformally invariant, so that components with this index position remain finite as r — > 00. 
This type of singularity is called "conformally compactifiable" [184] or "isotropic" [89]. For a 
classification of all possible global structures of spherically symmetric self-similar spacetimes see 
[103]. 

The global structure of the scalar field critical solution was determined accurately in [149] 
by assuming analyticity at the centre of spherical symmetry and at the past light cone of the 
singularity (the self-similarity horizon, or SSH). The critical solution is then analytic up to the 
future lightcone of the singularity (the Cauchy horizon, or CH) . Global adapted coordinates x and 
T can be chosen so that the regular centre r = 0, the SSH and the CH are all lines of constant 
X, and surfaces of constant r are never tangent to x lines. (A global t is no longer a global time 
coordinate). This is illustrated in Fig.]!] 
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Approaching the Cauchy horizon, the scalar field oscillates an infinite number of times but with 
the amplitude of the oscillations decaying to zero. The scalar field in regular adapted coordinates 
{x, t) is of the form 

r) ~ F,cg(T) + |x|^Fsi„g [r + //(r) + Kin , (30) 

where Frog(T), Fsing(T) and H{t) are periodic with period A, and the SSH is at a; = 0. These 
functions have been computed numerically to high accuracy, together with the constants K and e. 
The scalar field itself is smooth with respect to t, and as e > 0, it is continuous but not differentiable 
with respect to x on the CH itself. The same is true for the metric and the curvature. Surprisingly, 
the ratio m/r oi the Hawking mass over the area radius on the Cauchy horizon is of order 10~^ 
but not zero (the value is known to eight significant figures). 

As the CH itself is regular with smooth null data except for the singular point at its base, it 
is not intuitively clear why the continuation is not unique. A partial explanation is given in [149] . 
where all DSS continuations are considered. Within a DSS ansatz, the solution just to the future 
of the CH has the same form as (|30p . ficg(T) is the same on both sides, but -Fsing(''') can be chosen 
freely on the future side of the CH. Within the restriction to DSS this function can be taken to 
parameterise the information that comes out of the naked singularity. 

There is precisely one choice of Fsing(T) on the future side that gives a regular centre to the 
future of the CH, with the exception of the naked singularity itself, which is then a point. This 
continuation was calculated numerically, and is almost but not quite Minkowski in the sense that 
m/r remains small everywhere to the future of the SSH. 

All other DSS continuations have a naked, timelike central curvature singularity with negative 
mass. More exotic continuations including further Cauchy horizons would be allowed kinematically 
[ID] but are not achieved dynamically if we assume that the continuation is DSS. The spacetime 
diagram of the generic DSS continued solution is given in Fig. |3| 



3.4 Near-critical spacetimes and naked singularities 

Choptuik's results have an obvious bearing on the issue of cosmic censorship. (See jl89l for a 
general review of cosmic censorship.) Roughly speaking, fine-tuning to the black hole threshold 
provides a set of data which is codimension one in the space of generic, smooth, asymptotically flat 
initial data, and whose time evolution contains at least the point singularity of the critical solution. 
The cosmic censorship hypothesis must therefore be formulated as "generic smooth initial data 
for reasonable matter do not form naked singularities" . Here we look at the relation between 
fine-tuning and naked singularities in more detail. 

Christodoulou [61j proves rigorously that naked singularity formation is not generic, but in a 
rather larger function space, functions of bounded variation, than one would naturally consider. In 
particular, the instability of the naked singularity found by Christodoulou is not differentiable on 
the past light cone. This is unnatural in the context of critical collapse, where the naked singularity 
can arise from generic (up to fine-tuning) smooth initial data, and the intersection of the past light 
cone of the singularity with the initial data surface is as smooth as the initial data elsewhere. It 
is therefore not clear how this theorem relates to the numerical and analytical results strongly 
indicating that naked singularities are codimension- 1 generic within the space of smooth initial 
data. 

First, consider the exact critical solution. The lapse a defined by (pO)) is bounded above and 
below in the critical solution. Therefore the redshift measured between constant r observers located 
at any two points on an outgoing radial null geodesic in the critical spacetime to the past of the 
CH is bounded above and below. Within the exact critical solution, a point with arbitrarily high 
curvature can therefore be observed from a point with arbitrarily low curvature. Next consider a 
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Figure 3: The spacetime diagram of all generic DSS continuations of the scalar field critical 
solution, from [149] . The naked singularity is timelike, central, strong, and has negative mass. 
There is also a unique continuation where the singularity is replaced by a regular centre except at 
the spacetime point at the base of the CH, which is still a strong curvature singularity. No other 
spacetime diagram is possible if the continuation is DSS. The lines with arrows are lines of constant 
adapted coordinate a;, with the arrow indicating the direction of d/dr towards larger curvature. 
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spacetime where the critical solution in a central region is smoothly matched to an asymptotically 
flat outer region such that the resulting asymptotically flat spacetime contains a part of the critical 
solution that includes the singularity and a part of the CH. In this spacetime, a point of arbitrarily 
large curvature can be seen from ^+with finite redshift. This is illustrated in Fig.[4l 

Now consider the evolution of asymptotically flat initial data that have been flne-tuned to the 
black hole threshold. The global structure of such spacetimes has been investigated numerically in 
|1071 174 [ 177 1 1171j . Empirically, these spacetimes can be approximated near the singularity by ([6]). 
Almost all perturbations decay as the singularity is approached and the approximation becomes 
better, until the one growing perturbation (which by the assumption of fine-tuning starts out small) 
becomes significant. A maximal value of the curvature is then reached which is still visible from 
^^and which scales as [8Tj 

Rim,^ CC [P - P*f^ ■ (31) 

Finally, consider the limit of perfect fine-tuning. The growing mode is then absent, and all other 
modes decay as the naked singularity is approached. Note that ([6]) suggests this is true regardless 
of the direction (future, past, or spacelike, depending on the value of x) in which the singularity 
is approached. This seems to be in conflict with causality: the issue is sensitively connected to 
the completeness of the modes Zi and to the stability of the Cauchy horizon and requires more 
investigation. 

A complication in the supercritical case has been pointed out in [171] . In the literature on 
supercritical evolutions, what is quoted as the black hole mass is in fact the mass of the first 
apparent horizon (AH) that appears in the time slicing used by the code (spacelike or null). The 
black hole mass can and generically will be larger than the AH mass when it is first measured 
because of matter falling in later, and the region of maximal curvature may well be inside the 
event horizon, and hidden from observers at (See Fig. [5] for an illustration). The true black 

hole mass can only be measured at , where it is defined to be the limit of the Bondi mass me 
as the Bondi time ub — >■ oo. This was implemented in |171j . Only one family of initial data was 
investigated, but in this family it was found that tob converges to 10^^ of the initial Bondi mass 
in the fine tuning limit. More numerical evidence would be helpful, but the result is plausible. As 
the underlying physics is perfectly scale-invariant in the massless scalar field model, the minimum 
mass must be determined by the family of initial data through the infall of matter into the black 
hole. Simulations of critical collapse of a perfect fiuid in a cosmological context show a similar 
lower bound (114j due to matter falling back after shock formation, but this may not be true for 
all initial data |151j . 

3.5 Electric charge 

Given the scaling power law for the black hole mass in critical collapse, one would like to know 
what happens if one takes a generic 1-parameter family of initial data with both electric charge 
and angular momentum (for suitable matter), and fine-tunes the parameter p to the black hole 
threshold. A simple model for charged matter is a complex scalar field coupled to electromagnetism 
with the substitution Va — > Va + ieAa, or scalar electrodynamics. (Note that in geometric units, 
black hole charge Q has dimension of length, but the charge parameter e has dimension 1/length.) 

Gundlach and Marti'n-Garcfa [102J have studied scalar massless electrodynamics in spherical 
symmetry perturbatively. Clearly, the real scalar field critical solution of Choptuik is a solution 
of this system too. In fact, it remains a critical solution within massless scalar electrodynamics in 
the sense that it still has only one growing perturbation mode within the enlarged solution space. 
Some of its perturbations carry electric charge, but as they are all decaying, electric charge is a 
subdominant effect. The charge of the black hole in the critical limit is dominated by the most 
slowly decaying of the charged modes. From this analysis, a universal power-law scaling of the 



17 



Figure 4: Conformal diagram of the critical solution matched to an asymptotically flat one. 
RC=regular centre, S=Singularity, SSH=self-similarity horizon. Curved lines are lines of constant 
coordinate t, while converging straight lines are lines of constant coordinate x. Let the initial 
data on the Cauchy surface CS be those for the exact critical solution out to the 2-sphcrc R, and 
let these data be smoothly extended to some data that arc asymptotically flat, so that the future 
null infinity ,y+cxists. To the past of the matching surface MS the solution coincides with the 
critical solution. The spacetime cannot be uniquely continued beyond the Cauchy horizon CH. 
The redshift from point A to point B is finite by self-similarity, and the redshift from B to C is 
finite by asymptotic flatness. 
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Figure 5: The final event horizon of a black hole is only known when the infall of matter has 
stopped. Radiation at 1 collapses to form a small black hole which settles down, but later more 
radiation at 2 falls in to give rise to a larger final mass. Fine-tuning of a parameter may result in 
mi {p — p*)'^ , but the final mass m2 would be approximately independent of p. 

black hole charge 

Q-{p~P*)' (32) 

was predicted. The predicted value 5 ~ 0.88 of the critical exponent (in scalar electrodynamics) 
was subsequently verified in collapse simulations by Hod and Piran [121j and later again by Petryk 
|166j . (The mass scales with 7 ~ 0.37 as for the uncharged scalar field.) No other type of criticality 
can be found in the phase space of this system as dispersion and black holes are the only possible 
end states, though black holes with \Q\ < M can be formed 166'. 

General considerations similar to those in Sec. 12.61 led Gundlach and Martin-Garcfa to the 
general prediction that the two critical exponents are always related, for any matter model, by the 
inequality 

(5 > 27 (33) 

(with the equality holding if the critical solution is charged), so that black hole charge can always 
be treated perturbatively at the black hole threshold. This has not yet been verified in any other 
matter model. 

3.6 Self-interaction potential 

An example of the richer phenomenology in the presence of a scale in the field equations is the 
spherical massive scalar field with a potential m^c/)^ [33] coupled to gravity: In one region of 
phase space, with characteristic scales smaller than 1/m, the black hole threshold is dominated 
by the Choptuik solution and type II critical phenomena occur. In another it is dominated by 
metastable oscillating boson stars (whose mass is of order 1/m in geometric units) and type I 
critical phenomena occur. (For the real scalar field, the type I critical solution is an (unstable) 
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oscillating boson star |175j while for the complex scalar field it can be a static (unstable) boson 
star [ng.) 

When the scalar field with a potential is coupled to electromagnetism, type II criticality is 
still controlled by a solution which asymptotically resembles the uncharged Choptuik spacetime, 
but type I criticality is now controlled by charged boson stars |166j . There are indications that 
subcritical type I evolutions lead to slow, large amplitude oscillations of stable boson stars |1341 
I135| 1166] and not to dispersion to infinity, as had been conjectured in |115j . Another interesting 
extension is the study of the dynamics of a real scalar field with a symmetric double- well potential, 
in which the system displays type I criticality between the two possible vacua [123] . 

3.7 Nonspherical perturbations: stability and angular momentum 

Critical collapse is really relevant for cosmic censorship only if it is not restricted to spherical 
symmetry. Martin-Garcfa and Gundlach [146] have analysed all nonspherical perturbations of the 
scalar field critical solution by solving a linear eigenvalue problem with an ansatz of regularity at 
the centre and the SSH. They find that the only growing mode is the known spherical one, while all 
other spherical modes and all non-spherical modes decay. This strongly suggests that the critical 
solution is an attractor of codimension one not only in the space of spherically symmetric data but 
(modulo linearisation stability) of all data in a finite neighbourhood of spherical symmetry. 

More recently, Choptuik and collaborators [52] have carried out axisymmetric time evolutions 
for the massless scalar field using adaptive mesh refinement. They find that in the limit of fine- 
tuning generic axisymmetric initial data the spherically symmetric critical solution is approached 
at first but then deviates from spherical symmetry and eventually develops two centres, each of 
which approaches the critical solution and bifurcates again in a universal way. This suggests that 
the critical solution has non-spherical growing perturbation modes, possibly a single I = 2 even 
parity mode (in axisymmetry, only m = is allowed). There appears to be a conflict between the 
time evolution results J52' and the perturbative results T46' , which needs to be resolved by more 
work (see Sec. 15. 2p . 

Perturbing the scalar field around spherical symmetry, angular momentum comes in to second 
order in perturbation theory. All angular momentum perturbations were found to decay, and a 
critical exponent /i ~ 0.76 for the angular momentum was derived for the massless scalar field in 
[84] . This prediction has not yet been tested in nonlinear collapse simulations. 

4 More spherical symmetry 

The pioneering work of Choptuik on the spherical massless scalar field has been followed by a 
plethora of further investigations. These could be organised under many different criteria. We 
have chosen the following rough categories: 

• Systems in which the field equations, when reduced to spherical symmetry, form a single 
wave-like equation, typically with explicit r-dependence in its coefficients. This includes 
Yang-Mills fields, sigma models, vector and spinor fields, scalar fields in 2-1-1 or in 4-1-1 and 
more spacetime dimensions, and scalar fields in a semi-classical approximation to quantum 
gravity. 

• Perfect fluid matter, either in an asymptotically flat or a cosmological context. The linearised 
Euler equations are in fact wave-like, but the full non-linear equations admit shock heating 
and are therefore not even time-reversal symmetric. 

• Collisionless matter described by the Vlasov equation is a partial differential equation on 
particle phase space as well as spacetime. Therefore even in spherical symmetry, the matter 
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equation is a partial differential equation in four dimensions (rather than two). Intuitively 
speaking, there are infinitely more matter degrees of freedom than in the scalar field or in 
non-spherical vacuum gravity. 

• Spherically symmetric nonlinear wave equations on (3+1) Minkowski spacetime, and other 
nonlinear partial differential equations which show a transition between singularity formation 
and dispersal. 

Some of these examples were constructed because they may have intrinsic physical relevance 
(semiclassical gravity, primordial black holes), others as toy models for 3+1-dimensional gravity, 
and others mostly out of a purely mathematical interest. Table[l]gives an overview of these models. 

4.1 Matter obeying wave equations 
4.1.1 2d nonlinear a model 

We have already discussed in subsection 13 . 61 the effects of adding a potential to the evolution of the 
scalar field. An alternative generalization is a modification of the kinetic term in the Lagrangian, 
with the general form 

\g^-Gij{4,^)^^^'^,i,' (34) 

for an N-dimensional vector field ix) with / = l..iV, and where Gij{4>^) is a fixed nonlinear 
function acting as a metric on the so-called target space of the fields (j)^ . Such a system is called a 
nonlinear a model, harmonic map or wave map. The fields (j)^ and metric Gjj are dimensionless, 
and this allows the introduction of dimensionless parameters in the system, which cannot be 
asymptotically neglected using the arguments of subsection 12.61 (Compare with the potential 
V{(j)^), which has dimensions (length)^^, and hence requires dimensionful parameters). 

The case A'^ = 1 gives nothing new, and so Hirschmann and Eardly (HE from now on) studied 
the N = 2 case with a target manifold with constant curvature |_119. . proportional to a real 
dimensionless constant k. Using a single complex coordinate (f> the action of the system can be 
written as 

For K > this system is equivalent to the problem of a real massless scalar field coupled to Brans- 
Dicke (BD) gravity (with the BD coupling constant given by Sujbd = —12 + k~^). Licbling and 
Choptuik |143| have shown that there is a smooth transition in the BD system from DSS criticality 
for low K (the fiat target space case k = is equivalent to a self-gravitating complex massless scalar 
field, whose critical solution is the original DSS spacetime found by Choptuik) to CSS criticality 
for larger k (the case k = 1 is equivalent to the axion-dilaton system, and has been shown to 
display CSS criticality in |106j ). 

Generalizing their previous results for k = |118}I117] . HE constructed for each k a CSS solution 
based on the ansatz ^(r, x) = e'^^'^(f>{x) for the critical scalar field. Studying its perturbations HE 
concluded that this solution is critical for k > 0.0754, but has three unstable modes for k < 0.0754 
and even more for k < —0.28. Below 0.0754 a DSS solution takes over, as shown in the simulations 
of Liebling and Choptuik |143| . and HE conjectured that the transition is a Hopf bifurcation, 
such that the DSS cycle smoothly shrinks with growing k, collapsing onto the CSS solution at the 
transition and then disappearing with a finite value of the echoing period A. 

The close relation between the CSS and DSS critical solutions is also manifest in the construc- 
tion of their global structures. In particular, the results of [117] and [ST] for the CSS k = and 
K = 1 solutions respectively show that the Cauchy horizon of the singularity is almost but not 
quite flat, exactly as was the case with the Choptuik DSS spacetime (see subsection I3.3p . 



R 



(1-k|0|2)2 



(35) 
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4.1.2 Spherical Einstein-5f7(2) sigma model 



This is an = 3 sigma model, and it also displays a transition between CSS and DSS criticality, 
but this is a totally diflFerent type of transition, in particular showing a divergence in the echoing 
period A [HS] , 

In a reduction to spherical symmetry, the effective action is 



where r is the area radius, and the coupling constant 77 is dimensionless. It has been shown 
(numerically [25] and then analytically [26]) that for < 77 < 1/2, there is an infinite sequence of 
CSS solutions labelled by a nodal number n, and having n growing modes. (The case 7/ — 0, 
in which the sigma field decouples from gravity, will be revisited below.) The n-th solution is 
always regular in the past light cone of the singularity, but is regular up to the future light cone 
only for r? < 7y„ where t^q — 0.0688, ryi ~ 0.1518 and < r?n+i < 1/2. For larger couplings an 
apparent horizon develops and the solution cannot be smoothly continued. These results suggest 
that (pQ is a stable naked singularity for rj < rjQ, and 4>i acts as a critical solution between naked 
singularity formation and dispersal for rj < 770 and between black hole formation and dispersal for 
r]Q < rj < rji. The numerical experiments agree with this scenario in the range < 77 < 0.14. Other 
CSS solutions of this system are investigated in [53], and the possibility of chaos in [182] . 

Aichelburg and collaborators |126[ 1137) have shown that for 77 > 0.2 there is clear DSS type 
II criticality at the black hole threshold. The period A depends on 77, monotonically decreasing 
towards an asymptotic value for 77 — > cxi. Interesting new behavior occurs in the intermediate 
range 0.14 < 77 < 0.2 that lies between clear CSS and clear DSS. With decreasing rj the overall 
DSS includes episodes of approximate CSS |183j . of increasing length (measured in the log-scale 
time r). As ry — > 7;c — 0.170 from above the duration of the CSS epochs, and hence the overall 
DSS period A diverges. For 0.14 < rj < 0.17 time evolutions of initial data near the black hole 
threshold no longer show overal DSS, but they still show CSS episodes. Black hole mass scaling is 
unclear in this regime. 

It has been conjectured that this transition from CSS to DSS can be interpreted, in the language 
of the theory of dynamical systems, as the infinite-dimensional analogue of a three-dimensional 
Shil'nikov bifurcation [138] . High-precision numerics in [3] further supports this picture: for rj > rjc 
a codimension-1 CSS solution coexists in phase space with a codimension-1 DSS attractor such 
that the (one-dimensional) unstable manifold of the DSS solution lies on the stable manifold of 
the CSS solution. For 77 close to 77c the two solutions are close and the orbits around the DSS 
solution become slower because they spend more time in the neighbourhood of the CSS attractor. 
A linear stability analysis predicts a law A ~ — I- log(77 — 77c) -I- b for some constant b, where A is 
the Lyapunov exponent of the CSS solution. For rj = rjc both solutions touch and the DSS cycle 
dissapears. 

4.1.3 Einstein- Yang-Mills 

Choptuik, Chmaj and Bizoh [SD] have found both type I and type II critical collapse in the spherical 
Einstein- Yang-Mills system with SU(2) gauge potential, restricting to the purely magnetic case, 
in which the matter is described by a single real scalar field. The situation is very similar to 
that of the massive scalar field, and now the critical solutions are the well-known static 77 = 1 
Bartnik-Mckinnon solution [8] for type I and a DSS solution (later constructed in f93^) for type II. 
In both cases the black holes produced in the supercritical regime are Schwarzschild black holes 
with zero Yang-Mills field strength, but the final states (and the dynamics leading to them) can 
be distinguished by the value of the Yang-Mills final gauge potential at infinity, which can take 
two values, corresponding to two distinct vacuum states. 




(36) 
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Choptuik, Hirschmann and Marsa [54] have investigated the boundary in phase space between 
formation of those two types of black holes, using a code that can follow the time evolutions for 
long after the black hole has formed. This is a new "type III" phase transition whose critical 
solution is an unstable static black hole with Yang-Mills hair jT2lll88) . which collapses to a hairless 
Schwarzschild black hole with either vacuum state of the Yang-Mills field, depending on the sign of 
its one growing perturbation mode. This "colored" black hole is actually a member of a 1-parameter 
family parameterized by its apparent horizon radius and outside the horizon it approaches the 
corresponding BM solution. When the horizon radius approaches zero the three critical solutions 
meet at a "triple point" . What happens there deserves further investigation. 

Millward and Hirschmann jlSO, have further coupled a Higgs field to the Einstein- Yang-Mills 
system. New possible end states appear: regular static solutions, and stable hairy black holes 
(different from the colored black holes referred to above). Again there are type I or type II critical 
phenomena depending on the initial conditions. 

It is known that the spherical critical solutions within the magnetic ansatz become more un- 
stable when other components of the gauge field are taken into account, and so they will not be 
critical in the general case. 



4.1.4 Vacuum 4+1 

Bizoh, Chmaj and Schmidt [TH] have found a way of constructing asymptotically flat vacuum 
spacetimes in 4-|-l dimensions which are spherically symmetric while containing gravitational waves 
(Birkhoff's theorem does not hold in more than 3-f 1 dimensions). Recall that in (3-t-l dimensional) 
Bianchi IX cosmology the manifold is x where the is equipped with an SU{2) invariant 
(homogeneous but anisotropic) metric 

ds^ ^ Ljaf + Ljal + Llal (37) 

where the Ci are the SU (2) left-invariant 1-forms, and the Li are functions of time only. Similarly, 
the spacetimes of Bizon, Chmaj and Schmidt are of the product form x where the Li now 
depend only on r and t. This gives rise to nontrivial dynamics, including a threshold between 
dispersion and black hole formation. With the additional C/(l) symmetry Li = L2 there is only 
one dynamical degree of freedom. At the black hole threshold, type II critical phenomena are seen 
with A ~ 0.49 and 7 ~ 0.3289. In evolutions with the general ansatz where all Li are different, 
the U{1) symmetry is recovered dynamically |19j. Furthermore, the fact that each ?7(l)-symmetric 
solution exists in 3 copies gives rise to additional critical surfaces, and this lead to the prediction 
of 2-mode-unstable solutions, which was verified numerically, and the conjecture that there is a 
countable family of DSS solutions with n unstable modes jl9j . A similar ansatz can be made in 
other odd spacetime dimensions, and in 8-fl dimensions type II critical behaviour is again observed 

m. 



4.1.5 Scalar field collapse in 2+1 

Spacetime in 2+1 dimensions is flat everywhere where there is no matter, so that gravity is not 
acting at a distance in the usual way. There are no gravitational waves, and black holes can only 
be formed in the presence of a negative cosmological constant. (See [37] for a review.) 

Scalar field collapse in circular symmetry was investigated numerically by Pretorius and Chop- 
tuik |168] . and Husain and Olivier |128j . In a regime where the cosmological constant is small 
compared to spacetime curvature they find type II critical phenomena with a universal CSS criti- 
cal solution, and 7 = 1.20 ± 0.05 [168'. The value 7 ~ 0.81 [128] appears to be less accurate. 

Looking for the critical solution in closed form, Garfinkle [75] found a countable family of exact 
spherically symmetric CSS solutions for a massless scalar field with A = 0, but his results remain 
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inconclusive. The 9 = 4 solution appears to match the numerical evolutions inside the past light 
cone, but its past light cone is also an apparent horizon. The q = 4 solution has three growing 
modes although the top one would give 7 = 8/7~1.14if only the other two could be ruled out 
[83] . An attempt at this |120j seems unmotivated. At the same time, it is possible to embed the 
A = solutions into a family of A < ones [62l [63l HI] , which can be constructed along the lines 
of Sec. 12.61 so that Garfinkle's solution could be the leading term in an expansion in e~'^(— A)~^/^. 

4.1.6 Scalar field collapse in higher dimensions 

Critical collapse of a massless scalar field in spherical symmetry in 5+1 spacetime dimensions was 
investigated in [80]. Results are similar to 3+1 dimensions, with a DSS critical solution and mass 
scaling with 7 ~ 0.424. Birukou et al [TTJ 1127) have developed a code for arbitrary spacetime 
dimension. They confirm known results in 3+1 (7 ~ 0.36) and 5+1 (7 ~ 0.44) dimensions, and 
investigate 4+1 dimensions. Without a cosmological constant they find mass scaling with 7 ~ 0.41 
for one family of initial data and 7 ~ 0.52 for another. They see wiggles in the InM versus 
ln(p — p*) plot that indicate a DSS critical solution, but have not investigated the critical solution 
directly. With a negative cosmological constant and the second family, they find 7 = 0.49. [27] 
has made a more precise determination: 7 = 0.4131 ± 0.0001. This was motivated by an attempt 
to explain this exponent using an holographic duality between the strong coupling regime of 4+1 
gravity and the weak coupHng regime of 3+1 QCD [S], which had predicted 7 = 0.409552. 

Kol [133] relates a solution that is related to the Choptuik solution to a variant of the critical 
solution in the black-string black hole transition, and claims to obtain analytic estimates for 7 and 
A. This has motivated a numerical determination of 7 and A for the spherical massless scalar field 
in noninteger dimension up to 14 [ 178} 128] . 

4.1.7 Other systems obeying wave equations 

Choptuik, Hirschmann and Liebling [51] have presented perturbative indications that the static 
solutions found by van Putten [185' in the vacuum Brans-Dicke system are critical solutions. They 
have also performed full numerical simulations, but only starting from small deviations with respect 
to those solutions. 

Ventrella and Choptuik I186| have performed numerical simulations of collapse of a massless 
Dirac field in a special state: an incoherent sum of two independent left-handed zero-spin fields 
having opposite orbital angular momentum. This is prepared so that the total distribution of 
energy-momentum is spherically symmetric. The freedom in the system is then contained in a 
single complex scalar field obeying a modified linear wave equation in spherical symmetry. There 
are clear signs of CSS criticality in the metric variables, and the critical complex field exhibits a 
phase of the form e*""^ for a definite lo (the Hirschmann and Eardly ansatz for the complex scalar 
field critical solutions), which can be considered as a trivial form of DSS. 

Garfinkle, Mann and Vuille [86] have found coexistence of types I and II criticality in the 
spherical collapse of a massive vector field (the Proca system) , the scenario being almost identical 
to that of a massive scalar field. In the self-similar phase the collapse amplifies the longitudinal 
mode of the Proca field with respect to its transverse modes, which become negligible, and the 
critical solution is simply the gradient of the Choptuik DSS spacetime. 

Sarbach and Lehner [174] find type I critical behaviour in q + 3-dimensional spacetimes with 
C/(l) X SO{q + 1) symmetry in Einstein-Maxwell theory at the threshold between dispersion and 
formation of a black string. 
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Table 1: Critical collapse in spherical symmetry 
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4.2 Perfect fluid matter 
4.2.1 Spherical symmetry 

Evans and Coleman [69] performed the first simulations of critical collapse with a perfect fluid with 
equation of state p = kp (where p is the energy density and p the pressure) for fc = 1/3 (radiation) , 
and found a CSS critical solution with a mass-scaling critical exponent 7 « 0.36. Koike, Hara and 
Adachi jl31j constructed that critical solution and its linear perturbations from a CSS ansatz as 
an eigenvalue problem, computing the critical exponent to high precision. Independently, Maison 
[145] constructed the regular CSS solutions and their linear perturbations for a large number of 
values of k, showing for the first time that the critical exponents were model-dependent. As Ori 
and Piran before [161] 1162] . he claimed that there are no regular CSS solutions for k > 0.89, but 
Neilsen and Choptuik [1521 1153] have found CSS critical solutions for all values of k right up to 1 , 
both in collapse simulations and by making a CSS ansatz. The difficulty comes from a change in 
character of the sonic point, which becomes a nodal point for k > 0.89, rather than a focal point, 
making the ODE problem associated with the CSS ansatz much more difficult to solve. Harada 
jllO] has also found that the critical solution becomes unstable to a "kink" (discontinuous at the 
sonic point of the background solution) mode for k > 0.89, but because it is not smooth it does 
not seem to have any influence on the numerical simulations of collapse. On the other hand, the 
limit k —^ 1 leading to the stiff equation of state p — p is singular in that during evolution the fluid 
4-velocity can become spacelike and the density p negative. The stiff fluid equations of motion are 
in fact equivalent to the massless-scalar field, but the critical solutions can differ, dependending on 
how one deals with the issue of negative density |30] . Summarizing, it is possible to construct the 
Evans-Coleman CSS critical (codimension 1) solution for all values < fc < 1. This solution can be 
identified in the general classification of CSS perfect-fluid solutions as the unique spacetime that is 
analytic at the center and at the sound cone, is ingoing near the center, and outgoing everywhere 
else [SHI HU Ho]. There is even a Newtonian counterpart of the critical solution: the Hunter (a) 
solution |112] . 

p = kp is the only equation of state compatible with exact CSS (homothetic) solutions for 
perfect fluid collapse [36] and therefore we might think that other equations of state would not 
display critical phenomena, at least of type II. Neilsen and Choptuik [152] have given evidence 
that for the ideal gas equation of state (from now on EOS) p = kpoe (where po is the rest mass 
density and e is the internal energy per rest mass unit) the black hole threshold also contains a 
CSS attractor, and that it coincides with the CSS exact critical solution of the ultrarelativistic 
case with the same fc. This is interpreted a posteriori as a sign that the critical CSS solution is 
highly ultrarelativistic p — {1 + e)po ~ epo ^ po, and hence rest mass is irrelevant. Novak |158] 
has also shown in the case fc = 1, or even with a more general tabulated EOS, that type II critical 
phenomena can be found by velocity-induced perturbations of static TOV solutions. A thorough 
and much more precise analysis by Noble and Choptuik [156] 1157) of the possible collapse scenarios 
of the stiff fc = 1 ideal gas has confirmed this surprising result, and again the critical solution (and 
hence the critical exponent) is that of the ultrarelativistic limit problem. Parametrizing, as usual, 
the TOV solutions by the central density pc, they find that for low-density initial stars it is not 
possible to form a black hole by velocity- induced collapse; for intermediate initial values of pc, it is 
possible to induce type II criticality for large enough velocity perturbations; for large initial central 
densities they always get type I criticality, as we might have anticipated. 

Noble and Choptuik [156] have also investigated the evolution of a perfect fluid interacting 
with a massless scalar field indirectly through gravity. By tuning of the amplitude of the pulse it 
is possible to drive a fluid star to collapse. For massive stars type I criticality is found, in which 
the critical solution oscillates around a member of the unstable TOV branch. For less massive 
stars a large scalar amplitude is required to induce collapse, and the black hole threshold is always 
dominated by the scalar field DSS critical solution, with the fluid evolving passively. 
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4.2.2 Nonspherical perturbations 

Non-spherically symmetric perturbations around the spherical critical solution for the perfect fluid 
can be used to study angular momentum perturbatively. All nonspherical perturbations of the 
perfect fluid critical solution decay for equations of state p = kp with k in the range 1/9 < k < 
0.49 [97] ■ and so the spherically symmetric critical solution is stable under small deviations from 
spherical symmetry. Infinitesimal angular momentum is carried by the axial parity perturbations 
with angular dependence 1 = 1. From these two facts one can derive the angular momentum scaling 
law at the black hole threshold [95l [99] 

L^ip-p.r, (38) 

which should be valid in the range 1/9 < k < 0.49. The angular momentum exponent fi{k) is 
related to the mass exponent "/{k) by 

Mfc) = (2 + Ai(fc))7(fc). (39) 

where 

is the growth or decay rate of the dominant / = 1 axial perturbation mode. In particular for the 
value fc = 1/3, where 7 ~ 0.3558, ^i = (5/2)7 ~ 0.8895. 

Ori and Piran I162J have pointed out that there exists a CSS perfect fluid solution for < fc < 
0.036 generalizing the Larson-Penston solution of Newtonian fluid collapse, and which has a naked 
singularity for < fc < 0.0105. Harada and Maeda [1091 1112] have shown that this solution has 
no growing perturbative modes in spherical symmetry and hence a naked singularity becomes a 
global attractor of the evolution for the latter range of fc. This is also true in the limit fc = 0, 
which can be considered as the Newtonian limit [1131 1111] . Their result has been confirmed with 
very high precision numerics by Snajdr [177] . This seems to violate cosmic censorship, as generic 
spherical initial data would create a naked singularity. However, the exact result ()40l) holds for any 
regular CSS spherical perfect fluid solution, and so all such solutions with fc < 1/9 have at least one 
unstable nonspherical perturbation. Therefore the naked singularity is unstable to infinitesimal 
perturbations with angular momentum when one lifts the restriction to spherical symmetry. 

In the early universe, quantum fluctuations of the metric and matter can be important, for 
example providing the seeds of galaxy formation. Large enough fluctuations will collapse to form 
primordial black holes. As large quantum fluctuations are exponentially more unlikely than small 
ones, P{S) ^ e"* , where 6 is the density contrast of the fluctuation, one would expect the spectrum 
of primordial black holes to be sharply peaked at the minimal S that leads to black hole formation, 
giving rise to critical phenomena [154] . See also [SOI I194j . 

An approximation to primordial black hole formation is a spherically symmetric distribution of 
a radiation gas {p = p/3) with cosmological rather than asymptotically flat boundary conditions. 
In [1541 1155] type II critical phenomena were found, which would imply that the mass of primordial 
black holes formed are much smaller than the naively expected value of the mass contained within 
the Hubble horizon at the time of collapse. The boundary conditions and initial data were refined 
in [1141 II5I] . and a minimum black hole mass of ~ 10""* of the horizon mass was found, due to 
matter accreting onto the black hole after strong shock formation. 

4.3 Collisionless matter 

A cloud of collisionless particles can be described by the Vlasov equation, i.e., the Boltzmann 
equation without collision term. This matter model differs from field theories by having a much 
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larger number of matter degrees of freedom: The matter content is described by a statistical dis- 
tribution f{x^^,py) on the point particle phase space, instead of a finite number of fields 4){x^). 
When restricted to spherical symmetry, individual particles move tangentially as well as radially, 
and so individually have angular momentum, but the stress-energy tensor averages out to a spher- 
ically symmetric one, with zero total angular momentum. The distribution / is then a function 
f{r,t,p^,L'^) of radius, time, radial momentum and (conserved) angular momentum. 

Several numerical simulations of critical collapse of coUisionless matter in spherical symmetry 
have been published to date, and remarkably no type II scaling phenomena has been discovered. 
Indications of type I scaling have been found, but these do not quite fit the standard picture of 
critical collapse. Rein et al. [172] find that black hole formation turns on with a mass gap that 
is a large part of the ADM mass of the initial data, and this gap depends on the initial matter 
condition. No critical behavior of either type I or type II was observed. Olabarrieta and Choptuik 
[159] find evidence of a metastable static solution at the black hole threshold, with type I scaling 
of its life time as in Eq. (|13p. However, the critical exponent depends weakly on the family of 
initial data, ranging from 5.0 to 5.9, with a quoted uncertainty of 0.2. Furthermore, the matter 
distribution does not appear to be universal, while the metric seems to be universal up to an overall 
rescaling, so that there appears to be no universal critical solution. More precise computations by 
Stevenson and Choptuik [ 1791 , using finite volume HRSC methods, have confirmed the existence 
of static intermediate solutions and non-universal scaling with exponents ranging now from 5.27 
to 11.65. 

Martfn-Garci'a and Gundlach [147] have constructed a family of CSS spherically symmetric solu- 
tions for massless particles that is generic by function counting. There are infinitely many solutions 
with different matter configurations but the same stress-energy tensor and spacetime metric, due 
to the existence of an exact symmetry: two massless particles with energy-momentum in the 
solution can be replaced by one particle with 2p'^. A similar result holds for the perturbations. 
As the growth exponent A of a perturbation mode can be determined from the metric alone, this 
means that there are infinitely many perturbation modes with the same A. If there is one growing 
perturbative mode, there are infinitely many. Therefore a candidate critical solution (either static 
or CSS) cannot be isolated or have only one growing mode. This argument rules out the existence 
of both type I and type II critical phenomena (in their standard form, i.e., including universality) 
for massless particles in the complete system, but some partial form of criticality could still be 
found by restricting to sections of phase space in which that symmetry is broken, for example by 
prescribing a fixed form for the dependence of the distribution function / on angular momentum 
L, as those numerical simulations have done. 

A recent investigation of Andreasson and Rein [B; with massive particles has confirmed again the 
existence of a mass gap and the existence of metastable static solutions at the black hole threshold, 
though there is no estimation of the scaling of their life-times. More interestingly, they show that 
the subcritical regime can lead to either dispersion or an oscillating steady state depending on the 
binding energy of the system. They also conclude, based on perturbative arguments, that there 
cannot be an isolated universal critical solution. 

More numerical work is still required, but current evidence suggests that there are no type II 
critical phenomena, and that there is a continuum of critical solutions in type I critical phenomena 
and hence only limited universality. 

4.4 Criticality in singularity formation without gravitational collapse 

It is well known that the Yang-Mills field does not form singularities from smooth initial conditions 
in 3-1-1 dimensions [55] j but Bizoh and Tabor [53] have shown singularity formation in A+l (the 
critical dimension for this system from the point of view of energy scaling arguments) and 5-1-1 
(the first supercritical dimension). In 5+1 there is countable family Wn of CSS solutions with n 
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unstable modes, such that Wi acts as a critical solution separating singularity (Wq) formation from 
dispersal to infinity. In 4+1 there are no self-similar solutions and the formation of singularities 
seems to proceed through adiabatic shrinking of a static solution. 

Completely parallel results can be found for wave maps, for which the critical dimension is 
2+1. For the wave map from 3+1 Minkowski to the 3-sphere Bizoh [M] has shown that there 
is a countable family of regular (before the Cauchy horizon) CSS solutions labeled by a nodal 
number n > 0, such that each solution has n unstable modes. Simulations of collapse in spherical 
symmetry [211 1144] and in 3D [142] show that n = is a global attractor and the n = 1 solution 
is the critical solution. See also [65l [66] for computations of the largest perturbation-eigenvalues 
of Wq and Wi. Again, for the wave map from 2+1 Minkowski to the 2-sphere generic singularity 
formation proceeds through adiabatic shrinking of a static solution [22] . 

These results have led to the suggestion in [2^ that criticality (in the sense of the existence of 
a codimension-1 solution separating evolution towards qualitatively different end states) could be 
a generic and robust feature of evolutionary PDE systems in supercritical dimensions, and not an 
effect particular of gravity. 

Garfinkle and Isenberg [85] examine the threshold between the round endstate and pinching off 
in Ricci flow for a familiy of spherically symmetric geometries on S^. They have found intermediate 
approach to a special "javelin" geometry, but have not investigated whether this is universal. 

A scaling of the shape of the event horizon at the moment of merger in binary black hole 
mergers is noted in [42], but is really a kinematic effect. 

4.5 Analytic studies and toy models 
4.5.1 Exact solutions of Einstein Klein-Gordon 

A number of authors have explored the possibility of finding critical phenomena with CSS (rather 
than DSS) massless scalar critical solutions. The Roberts 1-parameter family of 3+1 solutions 
[173] has been analyzed along this line in [1631 [311 11911 1130] . This family contains black holes 
whose masses (with a suitable matching to an asymptotically flat solution) scale as {p — 1)-'^/^ for 
p > 1, but such a special family of solutions has no direct relevance for collapse from generic data. 
Its generalization to other dimensions has been considered in 72 . A fully analytic construction of 
all (spherical and nonspherical) linear perturbations of the Roberts solution by Frolov [7D1 [71] has 
shown that there is a continuum of unstable spherical modes filling a sector of the complex plane 
with ReX > 1, so that it cannot be a critical solution. Interestingly, all nonspherical perturbations 
decay. 

Frolov [72 has also suggested that the critical {p — 1) Roberts solution, which has an outgoing 
null singularity, plus its most rapidly growing (spherical) perturbation mode would evolve into the 
Choptuik solution, which would inherit the oscillation in t with a period 4.44 of that mode. 

A similar transition within a single 1-parameter family of solutions has been pointed out in 
pJO] for the Wyman solution [192] . 

Hayward [116[ [64] and Clement and Fabbri [62l [63] have also proposed critical solutions with 
a null singularity, and have attempted to construct black hole solutions from their linear pertur- 
bations. This is probably irrelevant to critical collapse, as the critical spacetime does not have an 
outgoing null singularity. Rather, the singularity is naked but first appears in a point. The future 
light cone of that point is not a null singularity but a Cauchy horizon with finite curvature. 

Other authors have attempted analytic approximations to the Choptuik solution. Pullin [170] 
has suggested describing critical collapse approximately as a perturbation of the Schwarzschild 
spacetime. Price and Pullin 169 have approximated the Choptuik solution by two flat space 
solutions of the scalar wave equation that are matched at a "transition edge" at constant self- 
similarity coordinate x. The nonlinearity of the gravitational field comes in through the matching 
procedure, and its details are claimed to provide an estimate of the echoing period A. 
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4.5.2 Toy models 

Birmingham and Sen |10j considered the formation of a black hole from the collision of two point 
particles of equal mass in 2+1 gravity. Peleg and Steif [165] have investigated the collapse of a dust 
ring. In both cases the mass of the black holes is a known function of the parameters of the initial 
condition, giving a "critical exponent" 1/2, but no underlying self-similar solution is involved. 

Frolov [75] and Frolov, Larsen and Christensen [75] consider a stationary 2+1-dimensional 
Nambu-Goto membrane held fixed at infinity in a stationary 3+1-dimensional black hole back- 
ground spacetime. The induced 2-t-l metric on the membrane can have wormhole, black hole, or 
Minkowski topology. The critical solution between Minkowski and black hole topology has 2-t-l 
CSS. The mass of the apparent horizon of induced black hole metrics scales with 7 = 2/3, superim- 
posed with a wiggle of period 37r/-\/7 in In p. The mass scaling is universal with respect to different 
background black hole metrics, as they can be approximated by Rindler space in the mass scaling 
limit. 

Horowitz and Hubeny |125j and Birmingham [9j have attempted to calculate the critical expo- 
nent in toy models from the adS-CFT correspondence. 

Burko [35] considers the transition between existence and non-existence of a null branch of 
the singularity inside a spherically symmetric charged black hole with massless scalar field matter 
thrown in. 

Wang [190] has constructed homothetic cylindrically symmetric solutions of 3-1-1 Einstein-Klein- 
Gordon and studied their cylindrically symmetric perturbations. It is not clear how these are 
related to a critical surface in phase space. 

4.6 Quantum effects 

Type II critical phenomena provide a relatively natural way of producing arbitrarily high curva- 
tures, where quantum gravity effects should become important, from generic initial data. Ap- 
proaching the Planck scale from above, one would expect to be able to write down a critical 
solution that is the classical critical solution asymptotically at large scales, as an expansion in 
inverse powers of the Planck length, see Sec l2.6l 

Black hole evolution in scmiclassical gravity has been investigated in 1+1 dimensional models 
which serve as toy models for spherical symmetry (see [88] for a review) . The black hole threshold 
in such models has been investigated in 43] , [7] , [181] , [130] , [195] , [1641 [25] . In some of these 
models, the critical exponent is 1/2 for kinematical reasons. 

In [M] a 3+1-dimensional but perturbative approach is taken. The quantum effects then give 
an additional unstable mode with A = 2. If this is larger than the positive Lyapunov exponent Aq, 
it will become the dominant perturbation for sufficiently good fine-tuning, and therefore sufficiently 
good fine-tuning will reveal a mass gap. 

5 Beyond spherical symmetry 

Numerical studies of critical collapse should go beyond spherical symmetry (and in the first instance 
to axisymmetry) for three reasons: 

• Weak gravitational waves in vacuum general relativity can focus and collapse. The black 
hole threshold in this process shows what in critical phenomena in gravitational collapse is 
intrinsic to gravity rather than the matter model. 

• Black holes are characterised by charge and angular momentum as well as mass. Angular 
momentum is the more interesting of the two because it is again independent of matter, but 
cannot be studied in spherical symmetry. 
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• Angular momentum resists collapse, but angular momentum in the initial data is needed to 
make a black hole with angular momentum. Therefore it is an interesting question to ask 
what happens to the dimensionless ratio J/Af^ at the black hole threshold. 

In the following we review what has been done so far. 



5.1 Perturbative approach to angular momentum 

We have already mentioned that when angular momentum is small, a critical exponent for |J| 
can be derived in perturbation theory. This has been done for the perfect fluid (see Sec. 14. 2p in 
first-order perturbation theory and for the massless scalar field (see Sec. 13. 7p . where second-order 
perturbation theory in the scalar field is necessary to obtain an angular momentum perturbation 
in the stress-energy tensor [84]. However, neither of the predicted angular momentum scaling laws 
has been verified in numerical evolutions. 

For a perfect fluid with equation of state p = kp with < A; < 1/9, precisely one mode that 
carries angular momentum is unstable, and this mode and the known spherical mode are the only 
two unstable modes of the spherical critical solution. (Note by comparison that from dimensional 
analysis one would not expect an uncharged critical solution to have a growing perturbation mode 
carrying charge.) The presence of two growing modes of the critical solution is expected to give rise 
to interesting phenomena jlOO) . Near the critical solution, the two growing modes compete. J and 
M of the final black hole are expected to depend on the distance to the black hole threshold and the 
angular momentum of the initial data through universal functions of one variable that are similar 
to "universal scaling functions" in statistical mechanics (see also the end of Sec. 12.71) . While they 
have not yet been computed, these functions can in principle be determined from time evolutions 
of a single 2-parameter family of initial data, and then determine J and M for all initial data 
near the black hole threshold and with small angular momentum. They would extend the simple 
power-law scalings of J and M into a region of initial data space with larger angular momentum. 



5.2 Axisymmetric vacuum gravity 

Abrahams and Evans [1] have numerically investigated black hole formation in axisymmetric vac- 
uum gravity. They write the metric as 



e^'^^^dr + dtf + r2e2''/3(d0 + dtf + g-^^/V^ sin^ 9 V] ^ (41) 



where the lapse a, shift components (3^ and and 3-metric coefficients and 77 are functions of 
r, t and 9. Axisymmetry limits gravitational waves to one polarisation out of two, so that there 
are as many physical degrees of freedom as in a single wave equation. On the initial slice, rj and Kg 
are given as free data, and i^, K'!^, and K'^ are determined by solving the Hamiltonian constraint 
and the two independent components of the momentum constraint. Afterwards, ry, Kg^ K"!^, and 
K"^ are evolved, and only is obtained by solving the Hamiltonian constraint, a is obtained by 
solving the maximal slicing condition Ki = for a, and and are obtained from the time 
derivatives quasi-isotropic spatial gauge conditions gee — r^grr and gre — 0. 

In order to keep their numerical grid as small as possible, Abrahams and Evans chose their 
initial data to be mostly ingoing. The two free functions rj and Kg in the initial data were chosen 
to have the same functional form they would have in a linearised gravitational wave with pure 
I = 2,m = angular dependence. This ansatz reduced the freedom in the initial data to one free 
function of advanced time. A specific peaked function was chosen, and only the overall amplitude 
was varied. 

Limited numerical resolution allowed Abrahams and Evans to find black holes with masses only 
down to 0.2 of the ADM mass. Even this far from criticality, they found power-law scaling of the 
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black hole mass, with a critical exponent 7 ~ 0.36. The black hole mass was determined from 
the apparent horizon surface area, and the frequencies of the lowest quasi-normal modes of the 
black hole. There was tentative evidence for scale echoing in the time evolution, with A ~ 0.6, 
with about three echos seen. Here rj has the echoing property ri{e^r, e^t) = r]{r,t), and the same 
echoing property is expected to hold also for a, 0, and r~^/3^. In a subsequent paper [2], some 
evidence for universality of the critical solution, echoing period and critical exponent was given 
in the evolution of a second family of initial data, one in which 77 = at the initial time. In this 
family, black hole masses down to 0.06 of the ADM mass were achieved. 

It is striking that at 14 years later, these results have not yet been independently verified. An 
attempt with a 3D numerical relativity code (but axisymmetric initial data), using free evolution 
in the BSSN formulation with maximal slicing and zero shift, to repeat the results of Abrahams 
and Evans was not successful [4]. The reason could be a combination of rather lower resolution 
than that of Abrahams and Evans, growing constraint violations, and an inappropriate choice of 
gauge. (It is now known that BSSN with maximal slicing and zero shift is an ill-posed system 

m-) 

5.3 Scalar field 

In 2003, Choptuik, Hirschmann, Licbling and Pretorius reported on numerical evolutions at the 
black hole threshold of an axisymmetric massless scalar field [52]. In axisymmetry with scalar 
field matter there is no angular momentum and only one polarisation of gravitational waves. The 
slicing condition is maximal slicing and the spatial gauge, in cylindrical coordinates, is gzz = dpp, 
Uzp — 0, similar to the gauge used by Abrahams and Evans. The Hamiltonian constraint is solved 
at every time step, and the time derivatives of the spatial gauge conditions are substituted into 
the momentum constraints to obtain second-order elliptic equations for the two shift components. 
Thus the evolution is partially constrained. Adaptive mesh-refinement was used in the numerical 
time evolution. 

The initial data were either time-symmetric or approximately ingoing, with the scalar field 
either symmetric or antisymmetric in z. In the symmetric case, even strongly non-spherical data 
were attracted to the known spherical critical solution for the massless scalar field. Scaling with the 
known 7 was observed in the Ricci scalar. However, with sufficiently good fine-tuning to the black 
hole threshold, the approximately spherical region that approaches the critical solution suffers an 
/ = 2 (and by ansatz m = 0) instability and splits into two new spherical regions which again 
approach the critical solution. The spatial separation of the two new centres is related to the 
smallest length scale that developed prior to the branching. There is evidence that with increasing 
fine-tuning each of these centres splits again. The antisymmetric initial data cannot approach a 
single spherical critical solution but the solution splits initially into two approximately spherical 
regions where the critical solution is approached (up to an overall sign in the scalar field). The 
separation of these initial two centres is determined by the initial data, but there is evidence that 
they in turn split. 

All this is consistent with the assumption that the spherical critical solution has, besides the 
known one spherical unstable mode, precisely one further I = 2 unstable mode. (Without the 
restriction to axisymmetry, if such a mode exists, it would be 5-fold degenerate with m — —2, . . . , 2.) 
This contradicts the calculation of the perturbation spectrum in [146) . Choptuik and co-workers 
do not state with certainty that the mode they see in numerical evolutions is a continuum mode, 
although they have no indication that it is a numerical artifact. The growth rate of the putative 
mode is measured to be A ~ 0.1 — 0.4, which should be compared with the growth rate A ~ 2.7 of 
the spherical mode and the relatively small decay rate of A ~ —0.02 claimed in |146j for the least 
damped mode, which is also an Z = 2 mode. We observe that the range of r in Figs. 6 and 7 of 
is about 10, and over this range the plot of the amplitude of the I = 2 perturbation against the 
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log-scale coordinate r seems equally consistent with linear growth in t as with exponential growth. 
(In the notation of [55], r denotes proper time and the accumulation point of echos, so that the 
log-scale coordinate r used in this review corresponds to — ln(r — t*) in the notation of [52 .) 

An interesting extension was made in [53] by considering a complex scalar field giving rise to 
an axisymmetric spacetime with angular momentum. The stress-energy tensor of a complex scalar 
field * is 

Tab - \ {"^M^fi + *,a*.b - .gab*,c*") (42) 

An axisymmetric spacetime with azimuthal Killing vector ^ admits a conserved vector field TabS,'', 
so that the total angular momentum 

J= f TabCn^dV (43) 

is independent of the Cauchy surface S. With the ansatz 

vE'(p,z,t,^) = <i>(p,z,i)e™^ (44) 

in adapted coordinates where ^ = d/dip and with m an integer, the stress-energy tensor becomes 
axisymmetric and hence compatible with the Einstein equations for an axisymmetric spacetime. 

For any E tangent to ^, in particular a hypersurface t = const., the angular momentum density 
measured by a normal observer becomes 

j = TabCn' = (n$ - n$) = mA^n'^S^a, (45) 

where $ = Ae^^ with S and A real and 11 = n°4>_a. By comparison the energy density measured 
by a normal observer is 

p = Tabu'^n^ = nn + Da^D''^, (46) 

where Da is the derivative operator projected into S. This means that the ratio of energy density 
to angular momentum density can be adjusted arbitrarily in the initial data, including zero angular 
momentum for a (f> that is real (up to a constant phase). On the other hand, even in the absence of 
angular momentum a purely real $ obeys a wave equation with an explicit m? / centrifugal term. 
For the same reason, regular solutions must have ~ on the axis, and there are no spherically 
symmetric solutions. Intuitively speaking, the centrifugal force resisting collapse appears unrelated 
to the angular momentum component of the stress-energy tensor in a way that differs from what 
one would expect in rotating fluid collapse or rotating (non-axisymmetric) vacuum collapse. 

For all initial data in numerical evolutions, a critical solution is approached that is discretely 
self-similar with log-scale period A ~ 0.42. (By ansatz this solution is axisymmetric but spherical 
symmetry is ruled out and so the critical solution cannot be the Choptuik solution.) The same 
critical solution is approached in particular for initial data with 11 = and hence no angular 
momentum, and initial data where 11 = $ and hence with large angular momentum. A scaling 
exponent of 7 ~ 0.11 is observed in the Ricci scalar in subcritical evolutions. The critical solution is 
purely real (up to an initial data-dependent constant phase) and hence has no angular momentum. 
Only m = 1 was investigated, but it is plausible that a different critical solution exists for each 
integer m. 

Far from the black hole threshold, J ^ AP in the final black hole, but nearer the black hole 
threshold, J ^ M^, where J and M are measured on the apparent horizon when it first forms. 
J/M"^ ^ is compatible with a non-rotating critical solution. 

In the absence of angular momentum, the wave equations for the real and imaginary part of <I> 
decouple. Assuming that the background critical solution is purely real with (5 = and A ~ 1, and 
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angular momentum is provided by a perturbation with 6 ~ e^'', one would expect p ~ | VAp ^ e^^ 
and j ~ A| VJ| ^ e'^^+'^''^. Integrating over a region of size e~^'^ when the black hole forms, we find 
M ~ e-^* and J ~ e(-2+^)^* . Then J ^ would imply A = -4. 

Lai [134] has studied type I critical phenomena for boson (massive complex scalar field) stars 
in axisymmetry, the first study of type I in axisymmetry. He finds that the subcritical endstate is 
a boson star with a large amplitude fundamental mode oscillation. 

5.4 Neutron star collision in axisymmetry 

A first investigation of type I critical collapse in an astrophysically motivated scenario was carried 
out in |129j . The matter is a perfect fluid with "Gamma law" equation of state P = (7 — l)pe, 
where F ~ 2 is a constant, P is the pressure, p the rest mass density and e the internal energy 
per rest mass (so that p{\ + e) is the total energy density). The initial data are constructed with 
P — kp^ for a constant k, which corresponds to the "cold" (constant entropy) limit of the Gamma 
law equation of state. The initial data correspond to two identical stars which have fallen from 
infinity. (The evolution starts at finite distance, with an initial velocity calculated in the first 
post-Newtonian approximation). The entire solution is axisymmetric with an additional reflection 
symmetry that maps one star to the other. 

The parameter of the initial data that is varied is the mass of the two stars. Supercritical data 
form a single black hole, while subcritical data form a single star. The diagnostics given are plots 
against time of the lapse and the Ricci scalar at the symmetry centre of the spacetime, and of 
outgoing I = 2 gravitational waves. Collapse of the lapse and blowup of the Ricci scalar are taken 
as indications of black hole formation. The limited numerical evidence is compatible with type I 
critical phenomena, with the putative critical solution showing oscillations with the same period 
in the lapse and Ricci scalar. For the critical solution to be exactly time-periodic, it would have 
to be spherical in order to not lose energy through gravitational waves, and there is some evidence 
that indeed it does not radiate gravitational waves. 

Other one-parameter families of initial data were obtained by fixing the mass of the stars in the 
first family very near its critical value and then varying one of the following: the initial separation, 
the initial speed, and F. For all approximately the same scaling law for the survival time of the 
critical solution was found. However, as all these initial data are very close together, this only 
confirms the validity of the general perturbation theory explanation of critical phenomena, but 
does not provide evidence of universality. 

A key question that has not been answered is how shock heating affects the standard critical 
collapse scenario. A priori the existence of a universal spherical critical solution is unlikely in the 
presence of shock heating as there is no dynamical mechanism to arrive at a universal spherical 
temperature profile, in contrast to non-dissipative matter models or vacuum gravity where all 
equations are wave-like. 

5.5 Black hole collisions 

Interesting numerical evidence for critical phenomena in the glancing collision of two black holes 
has been found by Pretorius [167] . He evolved initial data for two equal mass nonrotating black 
holes, with a reflection symmetry through the orbital plane, parameterised by their relative boost 
at constant impact parameter. The threshold in initial data space is between data which merge 
immmediately and those which do not (although they will merge later for initial data which are 
bound). The critical solution is a circular orbit that loses 1 — 1.5% of total energy per orbit through 
radiation. On both sides of the threshold, the number of orbits scales as 

n~ -7ln|]3-p*|, (47) 
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for 7 ~ 0.31 — 0.38. The simulations are currently limited by numerical accuracy to n < 5 and 
/S.E/E < 0.08 but Pretorius conjectures that the total energy loss and hence the number of orbits 
is limited only by the irreducible mass of the initial data, a much larger number. In particular, 
he speculates that for highly boosted initial data such that the total energy of the initial data is 
dominated by kinetic energy, almost all the energy can be converted into gravitational radiation. 

The standard dynamical systems picture of critical collapse, with the critical solution an at- 
tractor in the threshold hypersurface, appears to be consistent with these observation. Pretorius 
compares his data with the unstable circular geodesies in the spacetime of the hypothetical rotating 
black hole that would result if merger occurred promptly. These give orbital periods, and their 
linear perturbations give a critical exponent, in rough agreement with the numerical values for the 
full black hole collision. 

Pretorius does not comment on the nature of the critical solution, but because of the mass 
loss through gravitational radiation it cannot be strictly stationary. The mass loss would be 
compatible with self-similarity, with a helical homothetic vector field, but exact self-similarity 
would be compatible with the presence of black holes only in the infinite boost limit. Therefore 
the critical solution is likely to be more complicated, and can perhaps be written as an expansion 
with an exactly stationary or homothetic spacetime as the leading term. 

Pretorius also speculates that these phenomena generalise to generic initial data with unequal 
masses and black hole spins which are not aligned. This seems uncertain, given the claim by Levin 
[139] that the threshold of immediate merger is fractal if the spins are not aligned, and that the 
system is therefore chaotic. However, Levin's analysis is based on a 2nd order post-Newtonian 
approximation to general relativity, in which there is no radiation reaction, while the rapid energy 
loss observed here may suppress chaos. Nevertheless, the phase space is much bigger when the 
orbit is not confined to an orbital plane, and so the critical solution observed here may not be an 
attractor in the full critical surface. 
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